Eigenvalue curves of asymmetric tridiagonal random matrices 
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Abstract 

Random Schrodinger operators with imaginary vector potentials are studied in dimen- 
sion one. These operators are non-Hermitian and their spectra lie in the complex plane. We 
consider the eigenvalue problem on finite intervals of length n with periodic boundary con- 
ditions and describe the limit eigenvalue distribution when n — > oo. We prove that this limit 
distribution is supported by curves in the complex plane. We also obtain equations for these 
curves and for the corresponding eigenvalue density in terms of the Lyapunov exponent and 
the integrated density of states of a "reference" symmetric eigenvalue problem. In contrast 
to these results, the spectrum of the limit operator in l 2 (Z) is a two dimensional set which 
is not approximated by the spectra of the finite-interval operators. 

AMS subject classification: 65F15, 65F22, 15A18, 15A52, 60H25, 82B44 

Keywords: random matrix, Schrodinger operator, Lyapunov exponent, eigenvalue distribu- 
tion, complex eigenvalue. 



1 Introduction 



Consider an infinite asymmetric tridiagonal matrix J with real entries on the main diagonal 
and positive entries pk and r/% on the sub- and super-diagonal. Cut a square block of size 
n = 2m + 1 with the center at (0, 0) out of J n and impose the periodic boundary conditions in 
this block. The obtained matrix has the following form 



( Q—m f—m 

P—m+1 Q-m+1 r -m+l 



PO 



To 



P- 



\ 



Pm—1 Qm—1 fm— 1 
Pm Qm 
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We show nonzero entries of J n only. 

When n is large, the spectrum of J n cannot be obtained analytically, except for some special 
choices of pk, qk, and r k- However, it can be easily computed for "reasonable" values of n using 
any of the existing linear algebra software packages. If the pk, qk, and rk are chosen randomly 
the results of such computations are striking. For large values of n, the spectra of J n lie on 
smooth curves which change little from sample to sample. This was observed by Hatano and 
Nelson g, [fij. 
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Figure 1: Spectra of J n (n = 201) where (a) all non-zero entries are drawn from Uni[0, 1]; and 
(b) the sub-diagonal and diagonal entries are drawn from Uni[0, 1] and super-diagonal entries 
are drawn from Uni[^, l|]. 

Fig. [I] shows spectra of two matrices J n of dimension n = 201. One matrix has its non- 
zero entries drawn from Uni[0, 1]Q Its spectrum is shown on plot (a). Note that this matrix 
is only stochastically symmetric. For a typical sample from Uni[0, 1], J n ^ Jj. Nevertheless 
its spectrum is real. The other matrix has its diagonal and sub-diagonal entries drawn from 
Uni[0, 1] and super-diagonal entries drawn from Uni[J, li]. Its spectrum is shown on plot (b). 

Fig. |l] is in a sharp contrast to our next figure. We took the two matrices of Fig. || and 
subtracted \ from all their sub- and super-diagonal entries, including the corner ones. Fig. 
H shows spectra of the obtained matrices. Note that these spectra have a two-dimensional 
distribution. As will soon become clear, the eigenvalue curves in Fig. [I] are due to the sub- and 
super-diagonal entries having the same sign. 

The class of random matrices ( |1.1[ ) was introduced by Hatano and Nelson in 1996 |], [1C[] . 
Being motivated by statistical physics of magnetic flux lines and guided by the relevant physical 
setup, they considered random non-Hermitian Schrodinger operators H(g) = + ig) 2 + V 
and their discrete analog ues J n in a large box with periodic boundary conditions and discovered 
an interesting localization - derealization transition. Hatano and Nelson also argued that 
the eigenvalues corresponding to the localized states are real and those corresponding to the 
delocalized states are non-real. Since then there has been considerable interest to the spectra of 
J n and their multi-dimensional versions in the physics literature. 

1 Uni[a, b] denotes the uniform distribution on [a, 6]. 
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Figure 2: Spectra of J n (n = 201) where (a) the sub- and super-diagonal entries are drawn from 
Uni[— ^, ^] and the diagonal entries are drawn from Uni[0, 1]; and (b) the sub-diagonal entries 
are drawn from Uni[— ^, and the diagonal and super-diagonal are drawn from Uni[0, 1] 



Analytic descriptions of the spectrum of the random matrices J n , in the limit n — > oo, 
were obtained independently and almost simultaneously in [Q, |!| and in our paper [p~5|] . In this 
paper we present complete proofs of the results stated in |15|] . Our results do not rely on a 
particular asymptotic regime (as in ||) or a particular probability distribution (see ||]). We 
explain the phenomenon of eigenvalue curves depicted in Fig. |l[ Under general assumptions 
on the probability distribution of the p^, q^, and r^, we prove that the eigenvalues of J n tend 
to certain non-random curves in the complex plane when n — > oo. We also obtain an analytic 
description for these curves and for the corresponding limiting distribution of eigenvalues. This 
limiting eigenvalue distribution may undergo a transition from a distribution on the real line, 
as in Fig. 0(a), to a distribution in the complex plane, as in Fig. ||(b), when parameter values 
of the probability laws of the matrix entries vary. 

Our results resemble those obtained in the 1960's for truncated asymmetric Toeplitz matrices 
[19, 11]. The eigenvalues of finite blocks of an infinite Toeplitz matrix are distributed along curves 
in the complex plane in the limit when the block size goes to infinity, see [23] for a survey. Of 
course, this resemblance is only formal. The class of matrices we consider is very different from 
Toeplitz matrices. 

It is interesting that the spectra of the finite matrices J n with random entries, even in the 
limit n — > oo, are entirely different from the spectrum of the corresponding infinite random 
matrix J = tridiag(pj, qj, rj) considered as an operator acting on l 2 (Z). Indeed, suppose that 
the non-zero entries of J are independently drawn from a three-dimensional distribution with a 
bounded support S. Then, the probability is one that for any given (p, q,r) £ScK 3 and for 
any n G N and e > we can always find in J a block of size n such that for all j within this 
block \pj — p\ < e, \qj — q\ < e, and \rj — r\ < e. By the Weyl criterion, this implies that with 
probability one the spectrum of J = tridiag(pj, qj, Tj) contains the spectra of tridiag(p, q, r) for 
all (p, q, r) E S. (This argument is well known in the spectral theory of random operators, see 



3 



e.g. [l?]].) Since for every (p,q,r) the spectrum of tridiag(p, q, r) is the ellipse 

{z:z=pe it + q + r- it , t€ [0,2vr]}, 

the spectrum of the infinite tridiagonal random matrix J has a two-dimensional support with 
probability one if S is sumently rich. That is why the eigenvalue curves of J n are surprising. 
This discrepancy between the spectra of J n and J was mentioned in the preliminary account of 
our work [Jll| and was one of our motivations for studying the eigenvalue distribution for J n . 
Nothing of this kind happens for periodic sequences {(pj,qj,rj)}- If {Pk}, {%} and {r^} have 
a common period, then the spectrum of J lies on a curve in the complex plane [|16| coinciding 
with the eigenvalue curves of J n in the limit n — > oo. The two-dimensional spectra are specific 
to "sufficiently rich" random non-selfadjoint operators (see recent works J5], || which contain a 
much more detailed analysis of the spectral properties of infinite tridiagonal random matrices). 

One significant consequence of the above mentioned discrepancy between the spectra of J n 
and J is that the norm of the resolvent (J n — zl n )~ l may tend to infinity as n — > oo even 
when z is separated from the spectra of J n . This aspect of instability inherent in non-normal 



matrices |20| , 22] was thoroughly examined for asymmetric Toeplitz matrices in [18] and for 



random bidiagonal matrices in [21]. We do not discuss it here. 



The rest of this paper is organized as follows. Our main results, Theorems 2.1 and |2.2|, a re 
stated in Section ||[ We also discuss there several corollaries of our theorems. Theorem |2.1| is 
proved in Section || and Theorem [2^ is proved in Section |I[ Our proofs use a number of results 
from the spectral theory of (selfadjoint) random operators. These are summarized in Appendix. 

Acknowledgements We thank E.B. Davies and L.N. Trefethen for very useful and inter- 
esting discussions, and M. Sodin for bringing our attention to paper [O]. 



2 Main results and corollaries 



To simplify the notation, we label entries of J n by (j, k) with j and k being integers between 1 



and n. We also set pk 



and rfc = —e Vk , so that J n takes the form 



( 



\ 



qi 



o'/l 



3& \ 



oVn-1 



o'ln 



o£n— 1 



(2.1) 



/ 



where only the non-zero entries of J n are shownP). The corresponding eigenvalue equation can 
be written as the second-order difference equation 



1 i>k-l ~ e^fc+i + qk4>k = zipk, l<k<n, 



(2.2) 



with the boundary conditions (b.c.) 



tpo = ip n and VnH 



1 



tpl- 



(2.3) 



2 Following the tradition, we have chosen the sub- and super-diagonal entries to be negative. Provided the sub- 
and super-diagonal entries are of the same sign, the choice of sign does not matter. 
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Our basic assumptions are:- 



Vk, Qk)}kLo is a stationary ergodic (with respect to translations k — > k + 1) 
sequence of 3-component random vectors defined on a common probability (2-4) 
space (CI, J 7 , P); 

E ln(l + \qo\), E^o ; Eryo ar e finite. (2-5) 

The symbol E stands for the integration with respect to the probability measure, E / = J Q fdP. 

We have chosen the off-diagonal entries of J n to be of the same sign. More generally, one 
could consider matrices whose off-diagonal entries satisfy the following condition: 

(*) the product of (k, k + 1) and (k + 1, k) entries is positive for all k = 1, 2, ... n — 1. 

Any real asymmetric "purely tridiagonal matrix satisfying (*) can be transformed into a sym- 
metric tridiagonal matrix. The recipe is well known: put ip k = w kVk i n ( |2.2j ) and choose the w k 
so that to make the resulting difference equation symmetric. This is always possible when (*) 
holds, and the weights w k are defined uniquely up to a multiple. 

The above transformation is the starting point of our analysis. We set wq = 1 and 



Wk 



(2.6) 



Eqs. (p.2j)-( pT3D are then transformed into 

-c fe _i93fc_i - c k (p k+1 + q k p k = z<p k , 1 < k < n, 



(2.7) 



where 



fn+1 



W\ 1 

tpi, = <P0, 

W n +1 W n 



Cfe 



(2.8 



(2.9) 



In the matrix form the transformation from ip to cp can be expressed as J n i— > W J n W 
with W = diag(tui, . . . w n ). Eqs. ( ]2.7|) - ( |2.8| ) can be rewritten as W~ l J n Wip = zip, where 
ip = (ipi, ... ,ip n ) T . The asymmetry in the transformed eigenvalue problem is due to boundary 
conditions Q2.8|). One can visualize this asymmetry: obviously, 



W- 1 J n W = H n + V n , 



(2.10) 



( 11 

-ci 

V 



where 



-ci 



-c n -i 



-C n -1 



and Vn 



-c w n = -c e 



In \ 



(2.11) 



(2.12) 
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H n is a real symmetric Jacobi matrix. Its eigenvalue equation is given by Eq. (2.7) with the 
Dirichlet b.c. 

VWi = ^o = 0. (2.13) 



V n , which is due to (2.8), is a real asymmetric matrix. If E (£q — %) 7^ then one of the 
two non-zero entries of V n increases, and the other decreases exponentially fast with n (with 
probability 1). 

Though one could deal directly with J n , we deal with H n + V n instead. We thus consider 
the asymmetric eigenvalue problem ( |2.7D - ( |2.8| ) as an exponentially large "perturbation" of 
the symmetric problem Q2.7| ), ( 2.13] ), This point of view allows us to use the whole bulk 
of information about the symmetric problem (2/7), (|2.13| ) in the context of the asymmetric 
problem ( |2.7| )-( |2~8| ), It is worth mentioning that the exponential rate of growth of V n is very 
essential: no interesting effects would be observed for sub-exponential rates. 



Following the standard transfer-matrix approach, we rewrite Q2.7|) as 

<fk+1 )=A k (^ ), whereA fc = lf Qk ~ Z j. (2.11) 

Vk J V <Pk-l J c k \c k 1 



and introduce the notation 



S n (z)=A n -...-A 1 . (2.15) 



In order to formulate our results, we need to recall the classical notion of the Lyapunov 
exponent associated with Eq. fl2.7|) : 

j(z) = lim -Elog||S„(*)[|. (2.16) 

It is well known that the limit in (|2.16j ) exists for every z € C and is non-negative. Obviously any 
matrix norm can be used in ( |2.16| ), as they all are equivalent. It is convenient for our purposes 
to use the following norm 

||M|| = max^ \M jk \. (2.17) 



k 

3 



Introduce 



and consider the curve 



g = ±E( m -to) (2.18) 



C = {z G C : 7(z) = \g\ }. (2.19) 

This curve separates the two domains 

D 1 = {z£C: j(z) > \g\ } and D 2 = {z G C : j(z) < \g\ } (2.20) 

in the complex plane. Note that D 2 may be empty for some values of g. In this case, £ is either 
empty as well or degenerates into a subset of M. In Section we prove the following 
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Theorem 2.1 Assume ( \2.^ ) - ft2.!\ ). Then for P- almost all u = , rjj i Qj)}j^=o the following 
two statements hold: 

(a) For every compact set K\ C Di\R there exists an integer number n±(Ki,u)) such that for 
all n > m(Kx,u) there are no eigenvalues of J n in K\. 

(b) For any compact set K2 C D2 there exists an integer number ri2(K2,Lu) such that for all 
n > U2(K2,uj) there are no eigenvalues of J n in K2. 



To proceed, we need to introduce another well studied function, the integrated density of 
states iV(A) associated with Eq. (|2.7|). Let 



N n (\) = — #{ eigenvalues of H n in (— 00; A) }. 

n 

Then 

JV(A) = lim N n {\). (2.21) 

n— >oo 

It is well known that under assumptions ( [2.4D - Q2.5D , the limit in Eq. ( 2.2 1| ) exists on a set of 



full probability measure and that -/V(A) is a non-random continuous function, see e.g. |g, 17]. 

It is a fact from spectral theory of random operators that j(z) and iV(A) are related via the 
Thouless formula || 17] 



/-too 
log \z - A] dN{\) - E logc . (2.22) 
-00 

According to this formula, j(z), up to the additive constant E logco = \ E (£0 + tjq), is the 
log-potential of dN(X): 

/+00 
log \z - A] dN(X). (2.23) 
-co 

Then £ is an equipotential line: 

£ = {z£C: $(z) = max(E^ , Er/o) }. (2.24) 

This equipotential line consists typically of closed contours £j. In turn, each contour consists of 
two symmetric arcs whose endpoints lie on the real axis. The arcs are symmetric with respect 
to the reflection zi->I. The domain D2 defined above is simply the interior of the contours £j. 

Part (b) of Theorem 2J. states that for almost all {(£?, rjj , qj^JLo the spectrum of J n is 
wiped out from the interior of every contour £j as n — > 00. Parts (a) and (b) together imply 
that for P-almost all {(Cj7VjyQj)}JLo ^ ne eigenvalues of J n in the limit n — > 00 are located on 
£UR. 

Our next result describes the limiting eigenvalue distribution on £ U M. Let dv Jn denote the 
measure in the complex plane assigning the mass 1/n to each of the n eigenvalues of J n . 
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Theorem 2.2 Assume ^2.\ ) - Then for P- almost all to = {(^j : Tjj i Qj)}j^=o the following 

statement holds: For every bounded continuous function f(z) 



where 



lim / f{z)du Jn {z)= / f(X)dN(X)+ / f(z(l))p(z(l))dl, 
JS Jc 



£ = {AeK: AeSuppaW, $(A + iO) > max(Ef , E770) }, 



(2.25) 



(2.26) 



m = 5F 

and is the arc-length element on C. 



+oo 



dN(X) 



X-z 



(2.27) 



This theorem is proved in Section f|. Of course the eigenvalue curve C and the density of 
eigenvalues p(z) on it can be found explicitly only in exceptional cases]]. However, one can 
infer from our theorems rather detailed general information about the spectra of J n in the limit 
n — > oo. 

To facilitate the discussion, let us replace our basic assumption (^J) by the following more 
restrictive but still quite general one: 



{(Qk, £fc, r /fc)}fcLo ^ s a sequence of independent identically distributed random 
vectors defined on a common probability space. 



(2.28) 



Under assumptions fl2,4| ) and (2.28), the Lyapunov exponent j(z) is continuous in z every- 
where in the complex plane, see e.g. M. Also, there exist positive constants Co and xq depending 
on the distribution of (Cj^VjiQj) such that for all \x\ > xq 



log |x| - C < 7(s) < log |x| + Cq. 



(2.29) 



These inequalities are obvious if the law F of distribution of (£,j,r]j,qj) has bounded support. 
If the support of F is unbounded, then ( |2.29| ) can be obtained using methods of [14]. The 
continuity of j(z) together with ( 2.29|) imply that: 

(a) £ is not empty if and only if min 1 . e /j7(x) < \g\; 

(b) C is confined to a finite disk of radius R depending on the distribution of (£j,r)j,qj). 

To describe C we notice that ^(x + iy) is a strictly monotone function of y > 0. This follows 
from the Thouless formula (2.22). Hence, if j(x + iy) = \g\ then 



< \g\ 



(2.30) 



and, vice versa, for each x satisfying ( 2.30| ) one can find only one non- negative y(x) such that 
z = x + iy(x) solves the equation 

7(*) = \9\- (2-31) 
3 One such case ^| is when £j = g and rjj = —g for all j, and the diagonal entries qj are Cauchy distributed. 
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Because of the continuity of 7(x), the set of x where ( |2,30| ) holds is a union of disjoint intervals 
\aj,a'j] with aj < a'-. Therefore £ is a union of disconnected contours Cj. Each Cj consists of 



two smooth arcs, Uj(x) and —yj(x), formed by the solutions of Eq. (2.31) when x is running over 
[a,j,a'j]. Apart from the specified contours, the set of solutions of Eq. ( 2.3l| ) may also contain 
real points. These are the points where j(x) = \g\. 

It is easy to construct examples with a prescribed finite number of contours. However we do 
not know any obvious reason for the number of contours to be finite for an arbitrary distribution 

of Vj, Qj)- 



According to ( |2.25| ), the limiting eigenvalue distribution may have two components. One, 
represented by the first term on the right-hand side in ( 2.25| ), is supported on the real axis. We 



call this component real. The other, represented by the second term, is supported by C. We 
call this component complex. 

The following statements are simple corollaries of our Theorems. Assume (^J) and ( |2.28| ) 
and consider stochastically symmetric matrices J n , i.e. E£j = Er/j. In this case g = and hence 
the curve C is empty. Therefore the limiting eigenvalue distribution has the real component 
only. This is surprising and does not seem to be obvious a priori. But even more surprising is 
that C remains empty for all 

Id < g^' = mmj(x). 

If the distribution of (£,-, rjj, qj) is such that the support of the marginal distribution of qj contains 
at least two different points then 7(2;) is strictly positive for all x £ R 17]. In this case, by 
the continuity of j(x) and ( |2.29|) , > 0. 
On the other hand, if 

M>s£ 2) = 7( x )- 

x'Gbupp dN 

then £ of ( |2.26| ) is empty and the limiting eigenvalue distribution has the complex component 
only. Obviously, if the law of distribution of (£,j,i]jiQj) nas unbounded support, then g^ = +00. 
If g£) < \g\ < g^\ then the real and complex components coexist. 

It is worth mentioning that the density p(z) of the non-real eigenvalues, see ( |2.27| ), is analytic 
everywhere on C except the (real) end-points of the arcs. (If the limit eigenvalue distribution 
has no real component than p{z) is analytic everywhere on C.) The behavior of p{z) near an 
end-point of an arc, aj say, depends on the regularity properties of N(\) at this point. If the 
density of states dN(X)/d\ of the reference equation (|2.7|) is smooth in a neighborhood of A = aj 
then p{z) has a finite limit as z approaches aj along the arc. If, in addition, aj belongs to both 
£ and then the tangent to the arc at a,- exists and is not vertical. In other words, if dN{\)/d\ 
is smooth in a neighborhood of a branching point A = aj the complex branches grow out of this 
point linearly. 



3 Eigenvalue curves 



In this section we prove Theorem 2T. According to ( |2.10| ) the eigenvalues of J n and H n + V n 
coincide. It is more convenient for us to deal with H n + V n and we thus consider the eigenvalue 



problem djHgJ). By (|2.14|) - (|2.15[) , we may write (tp k+1 ,(p k 



k 
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1, . . . , n, instead of (2.7). In particular 



S n (z) 



<+>0 



On the other hand, as required by (|2T 



Therefore the eigenvalue problem (| 



1 



liU 



— 7 - B n S n (z 

Wn 



|2.8| ) is equivalent to the following one: 
= 0. 



Here I is 2 x 2 identity matrix and 

5 n = diag(/? n , 1), /3 n = = e *0»-eo+*»-*.). 

Thus z is an eigenvalue of if n + if and only if 



(3.1) 



(3.2) 



where fj^\z), i = 1,2 are the eigenvalues of B n S n {z). 

Without loss of generality we may suppose that g > 0. Then, by the ergodic theorem, 

1 1 Prob.l 

lim - log — = g > 0. (3.3) 

n— >oo n U7 n 

Since det S n (z) = co/c n , det B n S n (z) does not depend on z and, because of ([T^)and (|2? 



lim - log det B n S n {z) = lim -(770 + %) 

n— >oo ?l ra— >oo 71 



0. 



Let r(B n S n (z)) be the spectral radius of B n S n (z), i.e. 

r(B n S n (z)) = max{|/j| : /j is an eigenvalue of B n S n (z) }. 



(3.4) 



(3.5) 



/,From (3^3) - ( |3.4| ) we deduce a necessary condition for z to be an eigenvalue of H n + V^. This 
condition applies to P-almost all uj = , i]j , q^j^Lo an d is as follows: If z is an eigenvalue of 
H n + V n and n > n(u>) then 



1 11 

-\ogr(B n S n {z)) = - log — 
n n w n 

It is clear that this condition is necessary but not sufficient. 



(3.6) 



We start with part (b) of Theorem 2A. Let K be an arbitrary compact subset of Z?2- We 
shall prove that for P-almost all uj there is an integer n (K,uj) such that for all n > n (K,io) 
equation ( |3.6D cannot be solved if z E K. 
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Note that r(B n S n (z)) < \\B n S n (z)\\ < ||5 n ||||5„(z)||. Therefore, 

-logr(B n S n (z)) < itog||B n || + ilog||S n (*)|| 
n n n 

= o(l) H — log ||5 n (z)||, for almost all u. 
n 



(3.7) 
(3.8) 



Note that we only have to consider the case when g > 0. For, if g = then D2 = and we 
have nothing to prove. Obviously, one can find e > such that 

sup 7(2) < g - e. 

z&K 



But then, by Theorem A. 4 in Appendix, we have, for P-almost all u), 



limsup < sup — log ||5 n (z)|| 1 < sup 7(2) < g — e 

n^oo {z£K n J zeK 



and hence, by (3.7) 



limsup I sup - logr(S n S' n (z)) }> < g - e. 



Thus, for P-almost all uj, there exists an integer n\{K, u) such that for all n > m(K,co) 

- logr{B n S n (z)) < g - J for all z£K. 
n 2 



(3.9) 



Relations (|3.9|) and ( |3^3| ) contradict equation (3^). This proves part (b) of Theorem 2A 



Now we shall prove part (a) of Theorem 2.1. For this we need the following general (deter- 
ministic) result. 



Lemma 3.1 Let q%, ■ ■ ■ ,q n be real and cq, . . . ,c n be positive. For any two complex numbers (j)Q 
and (pi define recursively a sequence of complex numbers 4>2, ■ ■ ■ , 4>n+i as follows: 



Ck4>k 



+ 1 



(q k - z)(f>k ~ Ck-i<f>k-x, k = 1,2, . . . ,n. 



(3.10) 



Denote by S n the 2x2 matrix which maps (4>i,(/)q) t into (</> n +i, <ft n ) T obtained as prescribed 
above, i.e. S n (z) = A n ■ . . . ■ A\, where is given by ( \2.lA ). Then for every z G C + the matrix 
S n has two linearly independent eigenvectors (u*, 1) T and (v*, 1) T such that 



Imu* < 



Im z 



\u < 



z \ _|_ C n~l 



c n Im z 



(3.11) 



and 



lmv*>0, \v*\< 



Im z 



(3.12) 
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Proof. Define u%, . . . , as follows: u\ = u and Uk+i = fk(uk), k = 1, . . . n, where 

q k - z c fc _i 1 



fk(u) 



Ck 



c k u 



Obviously, u n+1 = F n (u) where F n (u) = / n (/n-i(- • • /i(«)) • • • )• 

Since Im z and all c k are positive, each of the functions fk maps C_ into itself and 



Tin 7 

Im/ n (n)< , |/„(«)|< 



In %n\ _|_ Cn—l 1 
Cn. Cn. U 



Vu G C. 



Therefore, F n maps continuously the compact set Q defined by the inequalities (3.11) into itself. 
By the fixed point theorem, there exists u* in Q such that u* = F n (u*), i.e. if ii n +i = u* given 
that u\ = u*. 

Now set cpo = 1, ipx = u* and iterate ( 3.10 ) to obtain (p^, k = 2, . . . n + 1 using this initial 
1 

¥>fc+l (Qk - z ) Cfc-l tPk-l 



data. Since ^- = u* and 

¥0 



k = 1,2 



Ck 



Ck fk 



we have that = Ufc+i f° r all with the as above. In particular, <p ™ +1 = u n+ \ = u* and 
hence ip n +i = fn,u*- But then 



1 



(3.13) 



and (it*, 1) T is an eigenvector of S n (z) with u* satisfying (^.11|) . 

To construct the other eigenvector one can iterate recursion ( |3. 10j ) in the opposite direction. 
Namely, set <p n +i = v n+ \, <p n = 1 and use ( 3.1CQ backwards to obtain the remaining (fk- Similarly 
to what we have done before, write ( |3.10| ) as Vk = fk( v k+i)i k = n, . . .1, where now Vk = 
and 



fk(v) 



Cfc-l 



qk- z- c k v 



Then v\ = F n (v n +±), where F n is the composition of fi, f%, . . . , f n . Each of these functions is 
continuous in the closure of the upper half of the complex plane and maps this set into itself. 
Since 



l/i(«)l < 



Im z 



Vn G C+, 



F n maps the compact set Q defined by the inequalities in ( p. 12 ) into itself. By the fixed point 
theorem, there exists v* G Q such that F n (v*) = v*, i.e. (pi = v*(fo given that (p n +i = v* and 
(p n = 1. But then 



^0 



(3.14) 



and (v*, 1) T is an eigenvector of S n (z) with v* satisfying ( |3 . 1 2[ ) . 

It is apparent that the two constructed eigenvectors are linearly independent. 
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Lemma 3.2 Let S n (z) be as in Lemma 3.1 and B n = diag(/? n , 1) with (5 n > 0. Then for every 
z £ C+ the matrix B n S n has two linearly independent eigenvectors (u n ,l) T and (v n ,l) T such 
that 



\mu„ < 



(5 n \mz 



\Un\ < 



0n 



+ 



and 



Cn Im z 



lmv n > 0, \v n \ < 



Im z 



(3.15) 
(3.16) 



Proof. B n S n (z) = A n ■ A n -i ■ . . . ■ A%, where At, . . . , A n -i as before (see ( p. 14 ) - ( p,15| )) and 

An. = 



Cn Cn 

1 







Since n > 0, Lemma 3.1 applies. 



We now return to our eigenvalue problem (|2.7|) - ( |2.8D and to the matrices S n (z) and B n 
associated with this problem, i.e. now are given by (^S|) and /3 n by Q3.2j). Set 



1 1 



(3.17) 



where (u n , 1) and (v n , 1) are the eigenvectors of B n S n (z) obtained in Lemma 3.2 . 



Lemma 3.3 Assume (2.4) - (2.1). Then for P-almost all {£kiVk>Qk}kLo the following state- 
ment holds: For all z € C+ 



lim -logll^llHT- 1 !! =0. 

n— >oo n 



(3.18) 



The convergence in ( SAh) is uniform in z on every compact set in C_ 



Proof. For any stationary ergodic sequence of random variables X n with finite first moment 
lim-n^oo X n /n = with probability 1. By ( 2.17| ), 

\2 



1 



< -logllTJIHT- 1 !! < -lo 



1 



Ki + \v n \ + 2) 



n 



n 



U n - V r , 



To complete the proof, apply inequalities ( |3.15| ) - ( p. 16 ). 

Recall that r(B n S n (z)) is used to denote the spectral radius (|3 



Lemma 3.4 Assume (2.4) - (2.1). Then for P-almost all {£k:Vk>Qk}kLo the following state- 
ment holds: For all z G C+ 



lim -logr(B n S n (z)) = j(z), 

n— >oo n 



(3.19) 



where j(z) is the Lyapunov exponent ( 2.16] ). The convergence in ( 3.18 ) is uniform on every 
compact set in C+ . 
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Proof. It follows from Lemma 3.2 that B n S n (z) = T n A n T n , where A n is the diagonal matrix 
of eigenvalues of B n S n {z) corresponding to T n . Then 



I T n [ 1 1 1 T n 



< 



\B n S n \ 



< l|r„||||r n 



-ii 



r(B n S n (z)) and, by Lemma EO, for P-almost 



With our choice ( p. 17 ) of the matrix norm, 
all {£k,Vk,Qk}kLo, 

lim — log — n ^ n ( = o uniformly in z on compact subsets of C_ 
n-*con r(B n S n (z)) 

On the other hand, for P-almost all {Ck,Vk, Qk}kLo, 



1 \\B n S n (z) 
hm -log 

n~*oo n \\b n {z)\\ 
This follows from the obvious inequalities 



uniformly in z. 



\S n (z 



\B r , 



£r < \\B n S n {z)\\ < \\B n \\\\S n (z)\\. 



Now the statement of Lemma follows from Theorem A.3 of Appendix. 



With Lemma "iA in hand, we are a in a position to prove part (a) of Theorem 2A. Let K 
be a compact subset of Z?i\R. As 7(2) is continuous in K and j(z) > g there, one can find an 
e > such that 



min7(z) > g + e. 

^From this, by Lemma |3.4| , for almost all u = {^, 9fc}^l , there exists an integer ni(K,uj) 
such that for all n > Ui{K, to) 



-\ogr{B n S n {z))>g+ £ - Mz^K. 
n 2, 



(3.20) 



In view of (|3.3j) , (|3.20[) contradicts (|3.6|) . Theorem |2T| is proved 



4 Distribution of eigenvalues 



We need to introduce more notations and to recall few elementary facts from potential theory. 

Let M n be an n x n matrix. We denote by dv Mn the measure on C that assigns to each of 
the n eigenvalues of M n the mass ^. This measure describes the distribution of eigenvalues of 
M n in the complex plane in the following sense. For any rectangle K C C 

v{K-M n )= [ dv Mn 
Jk 

gives the proportion of the eigenvalues of M n that are in K. The eigenvalues are counted 
according to their multiplicity. 



14 



The measure dv Mn can be obtained from the characteristic polynomial of M n as follows. Let 

p(z;M n ) = -log|det(M n -zI n )| (4.1) 



log|z-C|d^„(C) (4-2) 



c 



In view of (4.2), p(z;M n ) is the potential of the eigenvalue distribution of M n . Obviously, 
p(z; M n ) is locally integrable in z. Then for any sufficiently smooth function f(z) with compact 
support 



Ic 



log\z-(\Af(z)d 2 z = lim / log \z - (\Af(z)d 2 z 

e|0 J\ z -c\>e 



C £ ^J\z-(\>e 

= 2tt/(C), 

by Green's formula. Hence 

±- / p(z;M n )Af(z)d 2 z = [ f{z)du Mn {z). (4.3) 

Here A is the two-dimensional Laplacian and d 2 z is the element of area in the complex plane. 
Both p(z; M n ) and dv Mn define distributions in the sense of the theory of distributions and Eq. 
(|4.3| ) can be also read as the equality du Mn (z) = ^Ap(z; M n ) where now A is the distributional 
Laplacian. More generally, it is proved in potential theory that, under appropriate conditions 
on du, dv(z) = i^Ap(z), where p{z) = f log \ z — (\dv Un (Q is the potential of dv. This Poisson's 
equation relates measures and their potentials. 

In this section we shall calculate the limit of du Jn for matrices (|2.1|) , obtaining the potential 
of the limiting measure in terms of the integrated density of states N(X) (see (|2.21|)). Our 



calculation makes use of relation ( 2.10 ) according to which the asymmetric matrix J n is, modulo 
a similarity transformation, a rank 2 perturbation of the symmetric matrix H n (see (|2.11| )). The 
low rank of the perturbation allows to obtain explicit formulas describing the change in location 
of the eigenvalues. In our case, 

det(J n - zl n ) = det(F n + V n - zl n ) = d(z; H n , V n ) det(i? n - zl n ), (4.4) 



were 



d(z;H n ,V n ) = detlln + VniHn-zIn)- 1 } (4.5) 
= (l + a n G„i)(l + 6 n Gi 

n) a rfin 

G u G nn (4.6) 

with a n and b n being the top-right and left-bottom corner entries of V n and Gi m standing for 
the (l,m) entry of [H n — z/ n ) _1 . 

One can easily obtain ( |4.6| ) from ( |4.5| ) with the help of a little trick. Write V n in the form 
V n = A T B, where A and B are the following 2 x n matrices: 

a n ... \ / ... 1 

... 1 J I b n ... 



Then the n x n determinant in fl4.5| ) reduces to a 2 x 2 determinant, as det[J n + A T B(H. n 
zln)^ 1 } = det[/2 + B(H n — zI n )~ 1 A T ] and the latter, being expanded, gives 
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Eqs. (4.4) - Q4.6| ) yield the following relationship between the potentials of dv Jn and dv H 



p(z; J n 



p{z; H n + V n ) 

p(z;H n ) + -log\d{z;H n ,V n )\. 
n 



(4.7) 



The measures dv Hn are all supported on the real axis where they converge to the measure dN(X) 
when n — > oo. This implies the convergence of their potentials to &{z) = J K log \z — X\dN(X) for 



all non-real z, see Theorem A2 in Appendix. &{z) is the potential of the limiting distribution 
of eigenvalues for H n . Thus the main part of our calculation of p(z;H n ) in the limit n — > oo 
is evaluating the contribution of ^ log \d(z; H n , V n )\ to p(z;J n ). We can do this for all non- 
real z lying off the curve C (see ( 2. 19| ) and ( 2.24| )). The corresponding result is central to our 
considerations and we state it as Theorem 4.1 below. Note that, by the Thouless formula ( 2.22j ), 
j(z) = \g\ is equivalent to &(z) = max(E£o, E770) and thus (cf. ( 2.20 )) 



D x = {z : $(2) > max(E^ , E?7o) } and D 2 = {z : $(z) < max( E£ , Ery ) }• (4.8) 



Theorem 4.1 Assume ^4)-^^). Then for P- almost all {£,j,rjj,qj}JL Q , 

lim p(z; J n ) = <$>(z) Vz e 2MR; (4.9) 

n— >oo 

Wm. p(z-J n ) = max(E£ , Et? ) Vz G D 2 \M. (4.10) 



The convergence in - ( 4-K\ ) is uniform in z on every compact set in D\\M and D-^ 

respectively. 



Proof. In view of Eq. ( |4.7| ) and Theorem A. 2 in Appendix, we only have to prove that, with 
probability one, for any compact sets K\ C Z?i\R and K2 C D 2 \R 

lim - log \d(z; H n , V n )\ = (uniformly in z G K x C Di\R), (4.11) 

n— >oo n 

lim - log \d(z; H„, V n )\ = max( E£ , E770) - (uniformly in z € K 2 C D 2 \R). (4.12) 

n— »oo ji 

Let us write d(z; H n , V n ) in the form 

d(z; H n , V n ) = a n G n i + b n G\ n + a n b n G n \G\ n + (1 - a n bnG\\G nn ) (4.13) 



and estimate the four terms on the right-hand side (r.h.s.) of ( |4.13 ). Recall that a n and b n are 
the corner entries of V n and the G's are the corner entries of (H n — zln)^ 1 . In particular, 

G ln = G nl = n"=i^ = e^^\ 
det(H n -zI n ) det(H n -zI n )' 



Under assumptions (2A) - (|2.5| ), on a set of full probability measure, 

a n = -e n ^ E K°-»»)+°(i)] , b n = -e~ n ^ E «o-»»)+°(i)] , (4.14) 
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^"j/fe+li) _ „»[£E(fi>-H»)+o(l)] 



when n — > oo and 
On the other hand, 

| det(# n - zl n )\ = e np ^ H ^ = e "[*W-K«;ffn)] ) 



where r(z;H n ) = p(z;H n ) — 3>(z). According to Theorem A. 2 in Appendix for P- almost all 
Rr'/., •'//!.? o 

lim r(z; H n ) = (uniformly in z on every compact set K C C\M). (4-15) 

n— +00 

Therefore, with probability one, 

\a n G nl \ = e "[ESo-*( Z )+ z(i)] and | 6nGlri | = e n[Er, -*( 2 )+o 2 (i)] ^ ^_ 16) 



where the o z (l) terms vanish when n — > oo uniformly in z on every compact set K C 

To estimate the third term in the r.h.s. of ( 4.13D , recall the Thouless formula ( |2.22| ). As the 
Lyapunov exponent 7(2) is non-negative everywhere in the complex plane and $(z), for every 
fixed Rez, is an increasing function of | Imz|, we have that 

*(*) > ~E(£ + r/ ) VzeC\lR. (4.17) 

Moreover, since $(z) is continuous in z off the real axis, 

M(K) :=mm[*(z)-±E(t o + m )]>0 

for any compact set X C C\K. From this and ( |4.16| ), 

\a n b n G ln G nl \ = e nimo+m)-2H*)+oz(i)] (418) 

< e -n[2M(K)+o z (l)] y zeK , (4.19) 

In other words, the third term on the r.h.s. in ( [4.13 ) vanishes exponentially fast (and uniformly 
in z on every compact set in C\R) in the limit n — > 00. 

The fourth term in ( |4.13j ) cannot grow exponentially fast with n. Nor it can vanish exponen- 
tially fast. Estimating it from above is simple. Since \Gjj\ < \ j^ z \ , j = and a n b n = e°^ n \ 
we have, for P-almost all qjYjLoi that 

o(n) 

\l - a n b n G u G nn \ < 1 + H VzeC\K (4.20) 

with the o(n) term being independent of z. 

Estimating |1 — a n b n GuG nn \ from below is less trivial. We do this with the help of the two 
Propositions stated below. The first one is elementary and the second one is a standard result 
from spectral theory of random operators |f|. 
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Proposition 4.2 Let C a = {z G C : a < argz < a + ir }, where < a < ir. Then 



min 1 — z\ = sin a 

zeC a 



Proposition 4.3 Under assumptions ( $-4 ) ~ (%-5)> there exists a set Qq C Q of full probability 
measure such that for every u = {qk,£,k, r lk}kLo £ 



lim [{H n - zl n )- 1 ] u = P 
n— too 11 J 



da(X; uj) 
X-z 



Vz G 



(4.21) 



where o~(X;uj), for every uj G VLq, * s a measure on IR with unit total mass. The convergence in 
is uniform on every compact set in C\R. 



Remark. For almost all realizations oj, the semi- infinite matrix 



qi -ci 
-ci <?2 



-c 2 



where Ck = e^ h+Tlk \ specifies a selfadjoint operator H + (oj) on Z 2 (Z + ). Proposition |4.3| is a 
consequence of the selfadjointness of H + (uj), and the measure cr(A;u;) is simply the (1,1) entry 
of the resolution of identity for H + (uj). 



Let us set 



a(z; H n 



arg [(H n - zl n ) l ] n , if lmz>0 
arg [{zl n - F n ) _1 ] n , iflmz<0. 



Gn = [(H n — zl n ) as a function of z, maps the upper (lower) half of the complex plane 

into itself, and so does G nn . Therefore, 



< a(z;H n ) < ir \/z G C\I 



and 



a(z; H n ) < arg(a n b n G n G nn ) < a(z; H n ) + ir Vz G 
Then by Proposition |4.2| 

|1 - a n b n GnG nn \ > sin a(z;H n ) Vz G 



Obviously, sina(z; H n ) = | ImGn|/|Gn| and by Proposition for every uj G Q,q and every 
compact set K C 



lim sina(z; H n ) 



Im.z\ do~{X;uj) 
\X-z\ 2 



da(X; uj) 
X-z 



(uniformly in z G K). 



(4.22) 
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The r.h.s. in (4.22) is positive and continuous in z E K. Therefore it is bounded away from zero 
uniformly in z E K. Hence, for every uo E VLq (set of full probability measure) and for every 
compact set K C C\M 

(4.23) 



liminf |1 — a n b n G\\G nn \ > C(K;uj) > (uniformly in z E K), 



where the constant C(K;lo) depends only on K and to. 

Now we are in a position to prove fllHD and flP^ ). Let z E K C L>i\R. Then 

> ££ and <&(z) > Et] Vz E K. 

Therefore, by ( 4.16| ), the first three terms on the r.h.s. in ( 4.1 3| ) vanish exponentially fast (and 
uniformly in z E K) in the limit n — > oo. But then, with probability one, 



limsup — log \d(z; H n , V n )\ < (uniformly in z E K), 

n— *oo Tl 



in view of ( 4.20 ), and 



liminf — log \d(z; H n , V n )\ > (uniformly in z E K), 

n— >oo n 

in view of ( 4.23 ). This proves ( [4.11| ). 

Let K be any compact set in I?2\M. Then, by the definition of D2 and ( 4.17| ), 

min(E£ , Er] ) < < max(E£o, Er/ ). 

At this point we may assume, without loss of generality, that E^o / Eryo- For, if E£o = Et]q 
then IS e m pty This follows from ( 4.171 ). 

Now, if E£ > Er/o then Et] < $(z) < E^ for all z E K and, by ( pip , flOf ) - Q ), 
the first term on the r.h.s. in ( [4.13| ) dominates the other terms. In this case, with probability 
one, 



lim -log \d(z;H n ,V n ) 

n—>oo n 



E^o — &( z ) (uniformly in z E K). 



Similarly, if E£ < E770, then it is the second term that dominates and, with probability one, 



lim -\og\d(z;H n ,V n ) 

n—foo n 



Et]q — &(z) (uniformly in z E K). 



Theorem [4 . 1| is proved. 



We shall now deduce from Theorem 4.1 the weak convergence of the eigenvalue distributions 
dv Jn to a limiting measure in the limit n — > 00. In doing this we shall follow Widom [24, ^] 
who proved that the almost everywhere convergence of the potentials p(z;M n ) ( [4.2j ) of atomic 
measures dv Mn implies the weak convergence of the measures themselves provided they are 
supported inside a bounded domain in the complex plane. Under assumptions (2.4) - ( |2.-5| ). 
the spectra of J n ( |2.1[) are not necessarily confined to a bounded domain. To extend Widom's 
argument to our case we estimate the contribution of the tails of dv Jn to the corresponding 
potentials. 
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Define 



p(z) =max[$(z), E£ , E % ]. (4.24) 

This function coincides with the r.h.s. of ( [4.9[) - fl4.10| ) and is continuous everywhere in the 
complex plane except may be the set E, see ( [2. 26] ). 

&(z) is a subharmonic function Q and so is p(z). For, the maximum of two subharmonic 
functions is subharmonic too. Therefore, Ap(z) is non-negative in the sense of distribution 
theory and ^Ap(z) defines a measure in C which we denote by dv, 

dv{z) = — Ap(z). (4.25) 

First, we prove that in the limit n — > oo the potentials p(z; J n ) converge to p(z), for P-almost 
all {Ck,Vk, Qk}kLo, in the sense of distribution theory. 

In the Lemma below Co(C) is the space of continuous on C functions with compact support, 
and Cq°(C) is the subspace of Co(C) of those functions which are infinitely differentiable in Rez 
and Imz. 



Lemma 4.4 Assume ( \2-4 ) - l\2.3{ ). Then on a set of full probability measure, for every f € 
Co(C) ; and in particular for every f £ C^°(C), 

lim / f(z)p(z;J n )d 2 z = [ f(z)p(z)d 2 z. (4.26) 



C 



Proof. Since p(z) is subharmonic, p(z) € L Ioc (C), so the integral on the r.h.s. of ( 4.26 ) makes 



sense. Let Lj = {z £ C : dist(z,M U C) < 8}. By Theorem 4T, on a set of full probability 
measure S7 , 

lim p(z; J n ) = p(z) uniformly in z on compact subsets of C\(M U £). 
n— >oo 

It follows from this that on the same set fio 

lim I f(z)p(z;J n )d 2 z= lim f f(z)p(z)d 2 z (4.27) 

for every continuous f(z) with compact support and for every 5 > 0. 
Since f(z) has compact support, f(z)p(z) € L 1 (C). Therefore, 

Ve > 35 > : f \f{z)p{z)\d 2 z < e (4.28) 

Jl s 

and to complete the proof we only need to show that the same is true for f(z)p(z; J n ), uniformly 
in n. More precisely, it will suffice to prove the following statement. On a set of full probability 
measure, for every continuous function / with compact support 

Ve > 35 > : lim sup / \f(z)p(z; J n )\d 2 z < e. (4.29) 
n-+oo J Lx 
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Obviously, ( p~27D t ogether with ( gjg ) - ( pL29D imply ([O 



To prove fl4.29 ), we break up p(z; J n ) into two parts: 

1 n 

p(z;J n ) = -} log\zj -z\ (4.30) 



n 

3=1 



- T2 l °g\ z i - z \ + - yZ ^og\zj - z\, (4.31) 



n — ' n 

\zj\<R \zj\>R 



where z\, . . . ,z n are the eigenvalues of J n and the summation in the two sums in ( [4.31 ) is over 



all eigenvalues of J n satisfying the inequalities \zj\ < R and \zj\ > R respectively. The first term 
in ( |4.31| ) is bounded from above but is unbounded from below due to the log-singularities at 
Zj. On the contrary, the second term is bounded from below, provided z is separated from the 
boundary of the disk \z\ < R, but may be unbounded from above due to large values of \zj — z\. 
We shall treat these two terms separately. 

The required estimate on the integral in ( |4.29| ) involving the first term in the break-up of 
p(z; J n ) ( 4.3l|) can be obtained using the property of local integrability of log \z\. Recall that a 



family of functions {h a (z)} a ^ is called uniformly integrable in z on a bounded set D C C if 
for every e > there exists a 5 > such that for every compact set S C D of area less than 5 



/ \h a (z)\d 2 z < e Va G A. 
Js 



It is a corollary of the local integrability of log | z \ that for every compact set K C C the family 
of functions {log \Q — z\}c_^k is uniformly integrable in z on bounded subsets of C. From this 
one immediately obtains 



Proposition 4.5 Let X-r(ICI) be the characteristic function of the disk |£| < R. For every 
compact set K C C and for every R > the family of functions 



1 n 

- T]xR(ICj!)iog \Cj - z\\ 

Jn>l,Ci,-, 



3=1 

is uniformly integrable in z on bounded subsets of C. 



It is now apparent that for every continuous / with compact support and for every R > 



Ve > 36 > : lim sup 



f(z)~ £ log 



\zj\<R 



d 2 z < e. 



(4.32) 



To obtain an appropriate upper bound on the second term in ( 4.31| ), note the following. If 
C is such that dist(£, spec J„) > 1 then log \z — C| > for every z G spec J n and 



0< Yl lo gl^'-Cl < log|det(J n -C/„)| 

\zj\>R 



< ^log( e ^+e^ + |^| + |C|). 

3=1 



(4.33) 
(4.34) 
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The latter inequality is due to the fact that for every matrix A = WAjkW"" k=1 , \detA\ < 



IL=i Ysk=i \Ajk\- 



By Theorem 2.1, all non-real eigenvalues of J n are in the vicinity of C for all sufficiently large 
n. £ is a smooth curve and any vertical line in the upper half of the complex plane intersects C 
only once there. Therefore the probability is one that, moving up along the imaginary axis say, 
we can find a Co = W such that dist(Co, spec J n ) > 1 for all n > uq. 

Note that 



E ^ 

\zj\>R 



\Zj - Z\ 



E log \ Z 3 ~ Col 



\Zj\>R 



\zj\>R 



log 



1 + 



Co - z 



Z j ~ CO 



Choose now rg so that r$ > 1 and the disk \z\ < tq contains both Co an d the support K of 
/. Set R = 4ro- Then for all \z\ < ro and for all Zj in the exterior of \z\ < R = 4ro we have 
ICo — z\ < 2ro, \zj — Col > 3ro and [zj — z\ > 3ro and hence 



- E 



\zj\>R 



< log 3 + 



- E lo S z i ~C( 

\zj\>R 

1 - 

< log 3 + - 1o § ( e5j_1 + e% + \qj\ + 



V |z[ < ro 
[by (H-O)]. 

By the ergodic theorem, under assumptions ( |2.4[ ) - ( |2.5[ ), on a set ill of full probability measure, 
1 n 

lim - y^log (e^- 1 + e v > + lad + r ) = E log (e ?0 + e Vl + \qi \ + r ) < +c 

n^oc n — ' 



-oo 



Therefore, on fii, 



lim sup < 

n— +oo 



— sup 

n \z\<r 



E ^ 

jzj|>-R 



< Const 



and, because the disk \z\ < r covers the support of /, 



Ve > 35 > : lim sup 



/(^ E 



\zj\>R 



d 2 z < e. 



(4.35) 



fl4.32j ) and (|4.35|) imply (|4.29| ). Lemma [O] is proved. 



Corollary 4.6 Assume (2J^j - (2Ji). Then with probability one, 



lim iiz/j n = dv 

n— >oo 



zn t/ie sense of weak convergence of measures. 
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Proof. Since the operation A is continuous on distributions, Lemma 4.4 implies that, on a set 
of full probability measure, 



du Jn( z ) = TT-ApO; J n ) -> — Ap(z) = du(z), as n -> oo, 
as distributions. To complete the proof, recall that a sequence of measures converging as distri- 



butions must be converging weakly [12]. 



It is apparent that Ap(z) = everywhere off a line consisting of two parts. One is the 
equipotential line C (|2.24 ) that separates the domains D\ and Di\ the other is E ( 2.26 ) which 
is made up of all points of Supp dN that do not belong to the interior of the closed contours of 
C. 

p{z) is continuous in the upper and lower parts of the complex plane. It also has continuous 
derivatives everywhere but on C and E. Its normal derivative has a jump when z moves from 
D2 to D\ in the direction perpendicular to C. It follows from this that the restriction dv c of 
du to C\R is supported on C and has there density p(z) with respect to the arc-length measure 
dl on C. The density equals the jump in the normal derivative of p(z) multiplied by A 
straightforward calculation gives 



1 
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dN(X) 



X 



On the other hand, the restriction dv R of dv to the real axis is supported on E and coincides 
there with dN. Therefore, for every bounded continuous function f(z) 



f(z)du(z) 



f{z)du H {z)+ / f(z)du u {z) 
Jc 

f(\)dN(X) + / f(z(l))p(z(l))dl. 



(4.36) 
(4.37) 



Taking into account that the weak convergence of measures dv n is equivalent to the convergence 
of J f(z)du n (z) on bounded continuous functions we obtain from ( 4.36j ) - ( 4.37 ) and Corollary 
the statement of Theorem 



A Appendix 

In our analysis of the non-self- adjoint eigenvalue problem ([2.7|)-( |2~g| ) we have used a number of 
results about the finite difference equation 

-c fc _i</?fc-i - c k <p k+ i + q k <p k = z<p k (A.l) 

with random coefficients c k and q^. For the sake of completeness, we reproduce here the corre- 



sponding formal statements. With the exception of Theorem A. 4 these results are well known in 
the theory of random selfadjoint operators. Their proofs together with references to the original 
publications can be found in books 17]. Some of the results are proved there under slightly 



less general assumptions than those used in this paper. However only minor adjustments are 
needed to extend the published proofs to the generality of our assumptions. 
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Given a sequence of 2-component vectors {(c^, qk)}j^L\, consider a sequence of tridiagonal 
symmetric matrices H n of dimension n, n = 1,2, . . . , denned on {(c&, Qfc)}^Li in the following 

way 



/ ft 



-ci 



\ 



-ci 



' • -C n -1 
"Cn— 1 <Zn 



Assume that 



oj = {(cj , qj)}j^i is a stationary ergodic (with respect to the translation k — > k + 1) ,^ . 
sequence of random vectors defined on a common probability space (Q,J-,P); 



Cfc > for all fe, and E ln(l + \qx\) < +oo, E | In ci | < +oo. 



(A.3) 



Here as before the symbol E denotes averaging ever the probability space. 

The matrix H n is symmetric and has real eigenvalues. Their empirical cumulative distribu- 
tion function is defined as 



N n (X,u) 



1 



n 



^{eigenvalues of H n in (— oo,A)}. 



Theorem A.l Assume ( A. I ) - &A.3[ ). Then there exists a continuous non-random function 
N(X) of real variable X such that for almost all sequences lu 

lim NJX,u) = N(X). 

In other words, on a set of full probability measure, the eigenvalue counting measures dN n (X,uj) 
converge weakly, as n — > oo, to the limiting measure dN(X). 

The limit distribution function, A(A), is called the integrated density of states (IDS) of Eq. 
Let 



p(z, H n 



log | A — z\ dN n (X,uj) 



(A.4) 



and 



log |A - z| dN(X). 



(A.5) 



Under assumptions ( A.2 ) - ( |A.3| ), the integral in ( |A.5| ) converges for every non-real z, hence 
&(z) is well defined off the real axis. On the real axis, the equality in (A.5) is understood in the 
following senseQ 



$(x) 



lim 

ej.0 



/+oo 
max[log \X-x\, -e _1 ] dN(X), 
-oo 



x G 



Under assumptions (|A,2| ) - ( A.3 ), the prelimit integral converges for every e > 
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with the convention that <3?(x) = — oo if the above limit is — oo[]. With this convection, $(z) is 
subharmonic in the complex plane. In particular, 3>(z) is upper semi-continuous. 



Theorem A. 2 Assume (A.&) - (A.S). Then the following is true for almost all to: For all 
non-real z 



lim p(z;H n ) = $(z). 

n—>oo 

The convergence in ( A.(\ ) is uniform in z on every compact subset o/C\I 
Let 

Jn(z,Uj) = ~log||S n (.z)||, 

n 

where S n (z) is as in ( |2.15| ). 

The Lyapunov exponent is defined as follows: 



(A.6) 



(A.7) 



7(2) = lim E7„(z,lj). 

n— +00 

The limit above exists for every z £ C. Any matrix norm can be used in flA.7| ), as they all are 
equivalent. Since ||5 n || 2 > | detS n | and detS n (z) = co/c n , we have that for every z G C, and in 
particular for every real z, 

7(2) > 0. 



Theorem A. 3 Assume ( A.3j) - (jA.$) . Then for almost all 00 the following is true: For every 
compact set K C C\R 

lim 7 n (z,tj) = 7(2;) uniformly in z G K. 



Remark. In contrast to non-real z, on the real axis 

lim j n (x,u;) = 7(2;, u), (A.8) 

The limit above exists for almost all pairs (x,ui) and Ej(x,uj) = 7(2). However, when x is fixed, 
the set fl x of those uj for which the limit in (|A.7|) exists depends on x, and P(n x€ ^ft x ) = (see 

0, 0)- 

For every cv, j n (z,uj) is a subharmonic function in the complex plane. j(z) is also subhar- 
monic in C. This property of the Lyapunov exponent is very useful, see H. We use it here to 



deduce the following corollary from Theorem A. 3 



'In fact, $(z) is finite everywhere in the complex plane, see (A..9) 
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Theorem A. 4 Assume fiA.ty - (jA.3[ ). Then for almost all u> the following is true: For every 
compact set K C C, and in particular for every compact set K C R, 

limsup < sup7 n (z,w) > < sup j(z). 

If, in addition, j(z) is continuous in K then 



limsup< sup[7 n (z,u;) - j(z)] > < 0. 



Remark. This theorem plays a crucial role in our proof of the fact that the eigenvalues of J n are 
wiped out, as n — > oo, from the interior of each contour of the curve C, see part (b) of Theorem 
2,1 . Actually, our proof of part (a) of Theorem 2A also applies to any compact subset of Z^W- 
However, this proof is based on Theorem A. 3 and cannot be applied to the whole interior of C 
as it contains intervals of real axis. 



Theorem A. 4 follows immediately from Theorem |A.3| and the following result from the theory 



of subharmonic functions: 



Theorem A. 5 [see ftTSjl , p. 150] Let uj ^ — oo be a sequence of subharmonic functions in C 
converging in the sense of distribution theory to the subharmonic function u. If K is a compact 
subset of C and f is continuous on K , then 



limsup < sup(uj — /) > < sup(u — /). 

j^oo I K ) K 



Theorem A. 6 (Thouless formula) For all z G C, 7(2) = &(z) — E logci. 



It is a corollary of the Thouless formula and the positivity of the Lyapunov exponent "f(z) 
that 

$(z) > E log c \/z G C. (A.9) 
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